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Polaron and bipolaron dispersion curves in one dimension for intermediate coupling 
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Bipolaron energies are calculated as a function of wave vector by a variational method of Gurari 
appropriate for weak or intermediate coupling strengths, for a model with electron-phonon inter- 
actions independent of phonon wave vectors and a short-ranged Coulomb repulsion. It is assumed 
that the bare electrons have a constant effective mass. A two-parameter trial function is taken for 
the relative motion of the two electrons in the bipolaron. Energies of bipolarons are compared with 
those of two single polarons as a function of wave vector for various parameter values. Results for 
effective masses at the zone center are also obtained. Comparison is made with data of other authors 
for bipolarons in the Hubbard-Holstein model, which differs mainly from the present model in that 
it has a tight-binding band structure for the bare electrons. 
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I. INTRODUCTION 



Many authors have studied 

l^g Ql.2 .3.4.5.6.7.8.9 .10.11.12.13.14.15 



Of 

small 



energies 
and 

1 6 .17.18.19.20.21.22.23.24.25.26.27.28.29.30.31.32.33.34.35.36.37 



38,39,40,41,42,43 bipolarons in various numbers of dimen- 
sions, and some have made calculations of bipolaron ef- 
fective masses at the band minimum i^^^^^i^^ii^ 7 . 
However, we are not aware of published calculations of 
bipolaron energies as a function of wave vector which 
extend to large wave vectors except in the case of 
small bipolarons when the electron-phonon coupling 
is strong i 26 i 27 ' 28 In this paper we take up a study of 
this problem in a one-dimensional model with local 
interactions by use of a variational method used by 
Gurari for the single-polaron problem4i The method is 
appropriate for intermediate electron-phonon coupling 
strengths. It was discussed in some detail in a review 
article by Frohlichj^ 5 . and, for large polarons, gives the 
same results for binding energies and effective masses 
as those obtained by Lee et a/j 46 > 47 using different 
approaches. The Hamiltonian we shall use for the 
bipolaron problem is similar but not identical to that 
of the Hubbard-Holstein model ) 22 ' 24 i 48 ' 49 and will be 
formulated in terms of center-of-mass and relative coor- 
dinates of two electrons rather than in terms of electron 
creation and annihilation operators. This permits us 
to follow the variational method for the single-polaron 
problem with only minor modifications. The biggest 
difference of our model from the Hubbard-Holstein 
model is our assumption of a constant bare mass for 
the electron. We have not yet found a way to apply the 
Gurari method to the Hubbard-Holstein model itself. 

The original motivation for this work was to help to 
find out whether large enhancements of electron-electron 
attractions mediated via phonons or other intermediate 
bosons predicted in the simplest perturbation approach 
to interaction s 50 i 51 i 52 i 53 i 54 may still occur when complica- 



tions due to intermediate coupling are included. For met- 
als with Fermi energies large compared with the energy 
of any boson mediating electron-electron attractions, 
Eliashberg theory has been used to show that the net 
effect of enhancements of attractions is unlikely to lead 
to high-temperature superconductivity at high currents 
either in three-dimensional^ or quasi one-dimensional 
systems.^ 3 However, for small Fermi energies, infinite en- 
hancements are predicted by the simplest perturbation 
approach at some drift velocities in one dimensioni 52 i 53 
Study of the bipolaron problem in one dimension should 
give some insight into how this conclusion is modified 
by effects beyond perturbation theory in the limit of low 
concentrations of electrons. Although we shall use the 
term "phonon" for the boson mediating the attraction, 
we have in mind that interactions mediated by plasmons 
may be of practical importance for understanding possi- 
ble high-temperature superconductivity associated with 
bipolarons. 

There have been many published studies of bipo- 
larons in one dimension (including studies of two-site 
models) ^Q- 18 - 19 - 22 - 24 - 25 - 26 - 27 ' 28 - 30 - 32 - 33 ' 34 - 35 - 36 - 37 - 39 - 41 - 42 

but these mostly concentrate on finding the energies 
for the ground state as a function of coupling strength 
and Coulomb repulsion. Some calculate the bipolaron 
effective mass at the bottom of the band, but most 
do not discuss how energies vary with center-of-mass 
wave vector well away from the band minimum, except 
for strong electron-phonon coupling. Hohenadler et 
al£^ give graphical results for the spectral function of 
bipolarons as a function of wave vector for various cases 
where the coupling is not very strong, and approximate 
E{k) curves can be deduced from these. An easier com- 
parison to make is with some unpublished calculations 
performed by S. El Shawish for the Holstein model for 
some parameters corresponding roughly to some of those 
we have used. These values correspond to high ratios of 
phonon energy to electronic transfer integral t. 
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In Sec. |TT] we introduce the Hamiltonian and our vari- 
ational method. Some numerical results are presented in 
Sec. IIII| and some discussion is given in Sec. HV] 



II. 



HAMILTONIAN AND VARIATIONAL 
METHOD 



With a notation somewhat similar to that of Ref. HI 
but modified to apply to one dimension and for short- 
range interactions, we write the Hamiltonian H for the 
bipolaron problem with constant bare effective masses in 
the form 

H = H x + H u + H p + H c _ p . (1) 

Here Hx is the center-of-mass kinetic energy given by 

Hx = -^ 2 X , (2) 

where the center-of-mass coordinate 

X = \{ Xl +x 2 ), (3) 

and x\, X2 are the coordinates of the two electrons; H u 
is the Hamiltonian for relative motion, with 



H u = -2V 2 U + W(u), 



where 



U = Xl — X2 

is the relative coordinate and 



Wiu) 



otherwise , 



(4) 
(5) 

(6) 



with a the lattice constant; the phonon Hamiltonian H p 
and the electron-phonon interaction H c _ p are given by 



(7) 



and 



iJ c _ p = -iV(a/L)i ^[2cos{\ku)e lkX a k ] + H.c, (8) 



where a k and a k are creation and destruction operators 
for phonons of wave number k and L is the length of the 
crystal. We use the usual reduced units, with units of en- 
ergy, length and mass equal to fku, (fi/2mw) 5 and 2m re- 
spectively, where u> is the phonon angular frequency and 
m is the bare electron mass. The form of the potential 
term in H u is similar to but not the same as that in the 
Hubbard model because, with the form used, two elec- 
trons within a unit cell do not always interact, but this is 
compensated by interactions between electrons which are 
in neighbouring cells but separated by less than a. We 



do not include any spin-dependent terms in the Hamilto- 
nian, and in the following we shall not include any terms 
involving electron-spin wave functions. For bipolarons we 
implicitly assume that the two electrons in the pair have 
opposite spin by choosing a wave function for relative 
motion which is even in the relative coordinates. 

For a given center-of-mass vector Q, we adopt a trial 
wave function ^ of the form 



* = L-^e^ x <t )Q {u)\{^{Q 1 k,u)x, 



(9) 



where x is the phonon or other boson vacuum, 

^(Q,k,u) = N{Q,k,u) (10) 
x 1 + L^c{Q,k)cos{\ku)e~ ikx a\ , 

where c(Q, k) are variational parameters, 4>q(u) is a nor- 
malised even function, 



\(j) Q (u)\ 2 du = 1, 



(11) 



depending on one or more parameters, and N(Q, k, u) is 
a normalisation constant given by 



N(Q,k,u) = [l + L" 1 |c(Q,fc)| 2 cos 2 (ifcM)]" 



1. 



. (12) 

With such a trial wave function, calculations similar 
to those in Sec. 4 of Frohlich's review article^ for single 
polarons give the expectation values of the different terms 
in H. After replacing sums by integrals with use of the 
replacement 



4^ 2tt 



7r I a 



dk, 



find 



where 



2- 



7r / a 



— 77 / a 



\c(Q,k)\ 2 h k dk, 



— tc I a 



\c/)q(u)\ 2 cos 2 {\ku) du : 



(13) 



(14) 



(15) 



tt / a 



7r / a 



*{Q,k) - c{Q,k)]h h dk, (16) 



{Hx) 



2tt 



x I a 



Qk\c(Q,k)\ 2 h k 



(17) 



— ix I a 



47T 



7X fa 



c(Q,k)\ 2 k 2 h k dk 



/ ' kk'\c(Q,k)\ 2 \c(Q,k')\ 2 h kkl dkdk' . 

OTT 7 J -ir /a 
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Here 



< kk' 



\(/)q(u)\ 2 cos 2 (\ku) cos 2 (±k'u) du (18) 



and 



(H u ) — E u 



4tt 



r/a 



— 7T / a 



c(Q,k)\ 2 k 2 (l-h k )dk, (19) 



where the first term is independent of the variational 
coefficients c(Q,k). Going back to sums over k in 
Eqs. (fig]), flfi]) , (fT7)) and by minimisation of (i?) 

with respect to c(Q, A:) and c*(Q, £;), we find 



where 



c(Q,fc) 



G k = -!- / 7 fcWKQ^')! 2 *' 



i-(g-/ife 1 G fe )fc + i^ 1 fe 2: 



(20) 



2V 2 a r /a 



k'hkk/ 



-*/« [1 - (Q - h- k }G h <)k> + i/i^fc' 2 
From Eqs. (15]) and (fig]) . 



For the purpose of determining the coefficients we will 
make the following approximation 



hkk.i — hi-hk' 



This gives 



c(Q,k) 



-2iVa2 



\-{Q-g)k + \h^k 2 ' 



where 
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7T / a 



k J-v/a [l-Sk' + ih^k' 2 ] 



2 dk 7 



(21) with 
rdk'. 



S = Q-g. 



(23) 



(24) 



(25) 



(26) 



From Eqs. (T]), (TTJI, $T<g$i, (jTTI) , HU) and (MJ), we find that 
the expectation value of H, which we write as E(Q), is 
hkk' = \hk + ^h k ' + \hk+y + \hk-k' - \- (22) given by 



J 



E(Q) = m 2 -9 2 )- 

2V i a 2 



2V 2 a 

7T 



— rr/a 1 Sk -\- 2^k ^ 

kk'(h kk > - 



dk 



h k h k i 



n/a J-tt/o, {I ~ Sk + ^^(l - Sk' + \h~Jk' 2 ) 2 



dkdk' + E u 



(27) 



Although the variational coefficients determined in this 
approximation do not represent the optimal choice, 
nonetheless the expectation value of H, E(Q), yields an 
upper bound to the exact energy. 

The part E u of the expectation value of H u in Eq. (jTSJ) 
which is independent of c(Q, k) can be written as 

E u = T u + V u . (28) 

Here the relative kinetic energy T u is given by 

1 f \f k \ 2 k 2 dk, (29) 



T„ 



2vr 



where f k is given by 

fk = 



(u)e lku du 



(30) 



In Eq. (I29p the integral over k is, in principle, unre- 
stricted, i.e. is an integral over k from — oo to oo. How- 
ever, we shall restrict our trial wave functions for relative 



electron motion to those that do not have Fourier compo- 
nents for k > k m = n/a, and for such trial wave functions 
we use an integral from — k m to k m . 
Writing 

/•oo 

\<j>(u)\ 2 e iku du = d k (31) 

> 

we have, for a trial function for which f k = for |fc| > k m , 



2tt 



fc m -i|fc| 



(km 



fy-iJk'+hkdk'. 



(32) 



Using 



-iku 



du = t— rsin(i|fc|a), 
\k\ 



(33) 
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from Eq. ((6|) we find that the potential term V u in 
Eq. |gH|) is given by 



2P 



2ir 



2A-„ 



V u = — d k [sm(±ka)/k]dk (34) 



2A-„ 



We restrict the integral to the limits shown in Eq 
because dk vanishes for |fc| > 2k m . Remembering that 



cos = i[cos(fcw) + 1], 

Eqs. (fT5|) and (|3"Tj) enable us to write 



(35) 



(36) 



/ifcfe/ is then determined from Eq. I]22p. 

We now consider a two-paramctcr trial function for the 
relative motion, with an assumed function in real space 
modified by replacement of all Fourier components with 
wave vectors of magnitude greater than k m replaced by 
zero. 



Pseudo real-space function for relative motion 

We consider a function which has Fourier transforms 
up to |fc| = k m of the same form, up to a proportionality 
factor, as the transforms of </>, where 



4>{u) = JV (i + kM) e - AM , 



(37) 



with b and A adjustable parameters, and No & normali- 
sation factor. If b in the trial function is small, then the 
maximum of cj> is at u = 0, whereas if b > A there are two 
maxima at finite For |fc| > k m , we assume that the 
Fourier transform of <fi is zero. Thus, using Eq. (|3"T)) . fk 
of Eq. (f3T)]) is given by 



fk 



2A , 2fe(A 2 -fc 2 ) 
(A 2 + fc 2 ) "I" (A 2 + fc 2 ) 2 



if \k\ < k rn , 
if |&| > kjn , 



where AT is 



normalisation factor. 



(1/2tt) |//c| 2 = 1, we find TV is given by 



N 2 



where 



N = (1/X + b/X 2 + b 2 /2X 3 y 



(38) 



From 



(39) 



(40) 



is the normalisation factor of a wave function given by 
Eq. (f37|) before imposing a restriction on the Fourier com- 
ponents for |fe| > k m , and 



B= iLj fok{2 > 



(41) 



where fok is the same as fk but with the normalisation 
factor replaced by No- The departure from unity of the 



integral of Eq. (|4"Tj) gives the fractional change of the 
square of the normalisation factor due to putting the 
Fourier transforms of <f> for |fc| > k m as zero. 

From Eqs. (fT5")) and (|3T|) . if one were to suppose that 
the assumption of zero Fourier transforms of <f> for |fe| > 
k m only affected hk for \k\ < k m via the change of nor- 
malisation factor, we would find 



N 2 



2\ 



26(4A 2 - k 2 ) 2b 2 (8X 3 - 6Xk 2 ) 



4X 2 + k 2 (4A 2 + fc 2 ) s 



(4A 2 + fc 2 ) 3 



(42) 

for | A; | < k m . We have verified that, if the integral in 
Eq. ((32)) is extended to be from -oo to oo then Eq. (|36|) 
gives the result of Eq. . 

The kinetic term T u in Eq. ((28j) is given by Eq. Ip9|l . 
If cf> were taken as in Eq. ([3T]) without removal of the 
Fourier components for |fc| > k m , then using Eqs. ([4]), ([6]) 
and (|37p . the potential term V u in Eq. (|2"5|) would be 



K = P 



\(f>\ 2 du 



PN 2 [(1/X + b/X 2 + b 2 /2X 3 )(l - er Xa ) 
-(ab/X + b 2 a/2X 2 - 6 2 a 2 /4A)e" 



-Aai 



(43) 



By comparison of numerical results obtained from 
Eq. (34]) and ([36]) with those from Eqs. ((HJ) and ((43|) . 
we can find the effect of putting Fourier components of 
cj) for |fc| > k m on hk and y u for given values of the pa- 
rameters. 



III. NUMERICAL RESULTS 

Marsiglio 25 - states that the physical region of the 
Hubbard-Holstein model requires a condition a 2 / K < U, 
where a is a factor multiplying local vibrational coordi- 
nates Xi in the electron-phonon interaction, and K ap- 
pears in the expression \Kx 2 for the vibrational poten- 
tial energy. If we identify Pa in our model with U, al- 
though this correspondence is not exact, the same type 
of condition would require that 



2V 2 < Pa. 



(44) 



However, in attractive Hubbard models (U < 0), the 
negative U is usually thought to be brought about by 
effects of electron-phonon interactions overcoming a pos- 
itive Hubbard U, and Marsiglio's condition would appear 
to imply that negative U Hubbard models are unphysi- 
cal. Therefore we shall not restrict our numerical work 
to the region defined by Eq. (I44p . although in cases in 
which we are using our model with short-range forces as 
a simple approximation for effects of long-range forces we 
shall need to take the condition more seriously. 

For our numerical calculations we consider two types 
of cases. First we choose parameters which may be appli- 
cable to oxidised atactic polypropylene (OAPP) ) 55 i 56 and 
make use of some of the parameters used previously in the 
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model for superconductivity in channels composed of ar- 
rays of quasi one-dimensional filaments^ As in Ref. l53l . 
we take the cross section of individual filaments to be 
0.25 nm 2 . However, because the periodic potential due 
to aligned dipoles near the strings of charges forming the 
filaments according to the model of Grigorov et aZ . 57 i 58 i 59 
for individual filaments may not have very deep minima, 
we assume here that the bare electron mass is m c and 
not 2m c as in Ref. [H. We consider two variants of pa- 
rameters which could be applicable to OAPP. One, as 
in Ref. l53l . where the excitations mediating the electron- 
electron attraction are plasmons, and the other where 
optical phonons of average energy 0.36 eV mediate the at- 
traction. For the plasmon-induced interaction, we should 
strictly have long-ranged forces, but we hope that the 
model used here with short-range forces will give a first 
approximation to the real situation. We choose 0.36 eV 
for a phonon energy because there are several branches of 
the phonon spectrum associated with C — H 2 and C — H 3 
stretching vibrations whose energies at long wavelengths 
lie between 0.35 and 0.37 eV^ 

For given values of V 2 , a, Q and P, we solve for g 
of Eq. (|2"5"j) and then minimise the total energy with re- 
spect to A and b, using our full expressions of Eqs. (|32"|) 
and (f3"6"]l for h k , Eq. (|3"^f| for V , and the normalisation 
factor of Eq. (j3"9")l for f k . For computational purposes, 
we made use of programs or modifications of programs 
from a book^ 

In the model for strings or nanofilaments of Refs. l57l . 
I581I591 . if there is a periodic potential acting on the elec- 
trons in the string, it is due to groups of about three 
aligned dipoles surrounding each electron, and so the pe- 
riod of the potential, or lattice constant, a, is equal to 
the inverse of the linear electron concentration c — nd 2 , 
where n is the three-dimensional concentration in the 
filament and d 2 is its cross section. Thus we may not 
be free to choose the carrier concentration and lattice 
constant independently. However, we note that a recent 
theoretical study of channels through films of oxidised 
atactic polypropylene making use of Bose condensation 
of bosons in an array of nanofilaments with an E{K) 
curve for bosons consisting of a combination of linear 
and quadratic terms^* 2 - as indicated to occur in studies 
of Cooper pair dispersion^ 3 did not assume a periodic 
potential acting on the electrons in the nanofilaments. 
Another constraint was imposed in Ref. 62 because it was 
thought that it was probably necessary for the Fermi en- 
ergy to be smaller than a quarter of the energy of the 
excitation mediating the electron-electron attraction in 
order to have the possibility of large enhancements of in- 
teractions at high drift velocities i 53 However, in view of 
the results of the present paper, these enhancements may 
not occur when calculations of electron-electron interac- 
tions go beyond second-order perturbation theory. 

Bearing in mind the possible constraints, and assuming 
d 2 = 0.25 nm 2 and m h = m c as discussed above, we find 
two values of a and the related carrier concentrations 
corresponding to our two possible choices of phonons 



or plasmons to mediate the attraction. In both cases 
we use parameters such that the ratio of a to the po- 
laron radius r p — (fi/2mw)5 satisfies a/r p = 4. Us- 
ing a value for the high-frequency dielectric constant of 
2.2^ we find, for plasmons mediating the attraction, 
that a/r p = 4 implies a — 0.53 nm, the linear concen- 
tration c = 1/a = 1.9 x 10 7 cm , the three-dimensional 
carrier concentration n within a filament is n = 7.6 x 10 21 
cm" 3 , the plasmon energy hw p calculated using a three- 
dimensional formula appropriate for not too long wave- 
lengths is fuv p = 2.2 eV, the polaron radius r p = 0.132 
nm, and, from a one-dimensional formula, the bare Fermi 
energy e F = 0.34 eV. The values of most of these quanti- 
ties are only slightly different from those used in Ref [H. 
For phonons of energy 0.36 eV, a/r = 4 implies a = 1.3 
nm, r p = 0.325 nm, n = 3.1 x 10 21 cm" 3 and e F = 0.056 
eV. 

We also consider a different type of case appropriate 
for a quantum wire of a crystalline material. If, e.g., 
we assume that m = 2m c and huj = 0.05 eV, then the 
polaron radius is 0.62 nm, and a ratio a/r p — 0.5 would 
then imply a plausible value of the lattice constant of 0.31 
nm. We note that, for such a small value of a/r p , there 
is no point in doing calculations for very large values of 
Q, since the type of variational method we are using will 
not be appropriate when the single-polaron energy lies 
more than hut above the bottom of the band4£ 

Note that, when the polaron or bipolaron energy above 
the bottom of the band becomes close to hui, the dis- 
crepancy in energy from the types of states more gener- 
ally discussed (see e.g. Refs. l26U27ll28ll37h becomes large. 
This is because our method requires one or two electrons 
(for polarons or bipolarons) whose average wave vector 
or centre-of-mass wave vector is equal to the centre- 
of-mass wave vector of the polaron or bipolaron to be 
present whatever the energy of the state, whereas more 
commonly used methods find the lowest energy of the 
electron-phonon system for a given centre-of-mass wave 
vector. Such wave functions have only a small elec- 
tron component at the wave vector concerned when the 
threshold energy for emission of phonons is approached. 
We think that our method is concerned with states of 
more physical interest than the states usually discussed 
in this energy region. These states, with the wave vec- 
tor provided by phonons, and electrons at the bottom of 
the band, do not help in describing what happens when 
a polaron or bipolaron is accelerated rapidly past the 
threshold for emission of phonons. 

Since we are using a/r p = 4 for both possible param- 
eter choices for OAPP, the same computer calculations 
can be used for both plasmon or phonon-induced interac- 
tions, with only the values of quantities obtained in real 
units being different. Figures [1] and [2] show values of the 
bipolaron energy E(Q), the energy 2E S of two widely 
separated polarons each with wave vector 0.5Q, and the 
parameters A and b in the bipolaron trial function, for 
two values of V 2 and two related values of P for each V. 



6 



as a function of Q/Q , where 



Q 



2k 



2n/a. 



(45) 



For small b, 1/A is the bipolaron radius in units of the 
polaron radius. 

Figures [3]^a) and (b) show results of similar calcula- 
tions as a function of Q for a single value of V 2 for a 
more limited range of Q for the much smaller value of 
a/r p — 0.5 which may be appropriate for a quantum 
wire of crystalline material. 

The figures also show the values of g of Eq. (|25[) and of 
the parameters A and b which are found in the numerical 
work. We notice that in all cases shown in Figs. [T|— [3l 
there is a monotonic rise of energies with Q. 

We have also considered two other types of trial func- 
tions for relative motion, a Gaussian type of function 
with a second parameter in the prefactor, and a wave 
function constant in /c-space for wave vectors with mag- 
nitudes between minimum and maximum values k l and 
k 2 . Both these types of trial function gave poorer results 
for the bipolaron energies than the function used here. 
Also, in many cases k-, turned out to be zero in the sec- 
ond type of function, and so our second parameter often 
did not improve matters. 



IV. DISCUSSION 

Bipolaron energies measured from the bottom of the 
bare bands for the Holstein model obtained from unpub- 
lished results of S. El Shawish for the case of fwj/t = 16, 
g 2 = 0.125 and hui/t — 4, g 2 — 0.5, corresponding to 
a — 4, V 2 — 0.125 and a = 2, V 2 = 0.5 in our model, are 
equal to -0.335hu; and -1.22fiw compared with -0.29hui 
and -0.92hu> obtained by us, and so our results are 13% 
and 25% smaller than those of El Shawish for the mag- 
nitudes of this quantity. The energies of two single po- 
larons in our model for the same parameters are -0.21fku 
and -0.64/ia;, and so our bipolaron binding energies are 
0.08ftw and 0.28hui. Average effective masses (strictly 
the reciprocals of the average reciprocal masses) of bipo- 
larons up to Q = 0.1Q in our model for the same pairs 
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FIG. 1: Bipolaron energies E(Q), the energy 2E S of two 
widely separated single polarons, the parameters A and b of 
the trial wave function for relative motion of Eq. (|37p . and the 
quantity g of Eq. (|25|l as a function of bipolaron wave vector 
Q, for V 2 = 0.125, a = 4, and (a) P = 0, (b) P = 0.125. Ener- 
gies are measured with respect to the bottom of the bare band, 
in units of the phonon energy. Q m is defined by Eq. (|45|l . 



with a single dispersionless phonon for values of V 2 a up 
to three (corresponding to a up to 1.5 in their notation). 
For V 2 a = 1 (a = 0.5) and V 2 a = 3 (a = 1.5), from 
Figs. 5 and 8 of their paper we find E w Q.52fku and 



of parameters are 2.13m b and 2.69m b , whereas the aver- E s w L69fru;, while their masses m s satisfy mjm h w 1.2 
age masses over the same region of wave vectors inferred 
from the results of El Shawish are 2.88m b and 3.05m b . 



Thus the departure of the mass from twice the single- 
particle mass is much larger in the results of El Shawish. 
For a — 4, with V 2 = 0.25 and 0.375, bipolaron binding 
energies in our model are Q.22hw and 0.38huj respectively, 
while bipolaron masses for these two cases are 2.31m b and 
2.50m b . Thus we are able to get bipolaron binding en- 
ergies of over a third of the phonon energy without very 
large effective masses. 

In the limit that a — > (with ta 2 constant), our method 
gives single-polaron binding energies E s = 0.5V 2 ahu) and 
effective masses m s such that m s /m h — 1 + 0.25V 2 a. 
These results may be compared with those of Chen et 
a/. 65 for a continuum model with short-range interactions 



and 3.7 for the two values of coupling. Thus E s for the 
two models differ by only 13% for V 2 a = 3, whereas 
masses for the two models are within 4% of each other 
for V 2 a = 1 (closer to 20% for the departure of m/m b 
from unity), but for V 2 a — 3 the discrepancy in masses is 
large (m s /m h — 3.7 in the model of Ref . Ifia. but 1.75 for 
our model). The larger discrepancy for masses between 
the two models than for binding energies is related to the 
fact that intermediate-coupling methods are accurate for 
binding energies up to higher values of coupling constants 
a in the Frohlich model than for masses i 47 i 66 

Another case for which we may make comparison with 
other authors is that of a = 1, V 2 = 1, corresponding to 
oj/t = 1 and g 2 = 1 in the Hubbard-Holstein model. Here 
we find that the bipolaron energy measured from the bot- 
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FIG. 2: As for Fig. [T] but with V 2 = 0.25, a = 4 and (a) 
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torn of the bare band is — 1.10/kj, and the energy of two 
single polarons is — 0.80/kj. Thus the bipolaron binding 
energy with respect to two single polarons is O.SOHlo in 
our model, compared with about 0.5/kj which can be in- 
ferred from the inset to Fig. 3 of a paper by Bonca et 
al 35 For this case we have also found the value of Pa 
in our model for which the bipolaron binding energies 
for the above values of parameters vanish. We find the 
bipolaron energy vanishes for Pa « 1.9. The result of 
a vanishing of the binding energy at Pa rj 1.9 may be 
compared with U ~ 1.6 for the binding energy to vanish 
which may be inferred from the inset to Fig. 3 of Ref. [35l . 

The smaller binding energies and masses, and larger 
partial bandwidths cf. what are probably accurate re- 
sults for the Hubbard-Holstein model are thought to be 
due to a combination of (i) limitations of our variational 
method to quite weak couplings, (ii) the use of only a 
two-parameter model for the relative motion, and (iii) to 
differences between a Holstein model and a model with 
constant bare mass. We do not know which of the three 
causes of discrepancy is most important at present. 

A lower limit for the percentage change due to depar- 
ture from weak coupling of the expectation value of the 
part of the Hamiltonian which does not depend on rel- 
ative motion of the two electrons can be obtained by 



looking at the percentage change in single-polaron ener- 
gies at the same value of coupling when this is known. 
An approximate upper limit to the percentage change 
in the same quantity can be obtained from the percent- 
age change in the single-polaron binding energy when 
the coupling g 2 is twice as large, since, in the limit when 
the relative wave function is very small, the coupling to 
phonons is twice as large for a bipolaron as for a polaron. 
For example, for the case g = 1, t/huj = 1, in the Hol- 
stein model, corresponding to V 2 = 1 and a = 1 in our 
model, one can estimate from Fig. 1 of Ref. [67| that both 
the lower and upper limits in the errors in the bipolaron 
energy below the bottom of the bare band due to use of 
a weak-coupling method are very small. 

A simple way to estimate errors due to the change 
from the constant bare-mass model to the Holstein model 
is to compare the weak-coupling result for the single- 
polaron model in the Holstein model given in Eq. (8) 
of Ref. 67J, which gives a polaron binding energy E hs of 
g 2 huj/(l + 4J/hu>)^ in their notation. 68 In our notation, 
J/hu corresponds to 1/a 2 , and g to V. Thus their result 
in our notation is 

E hs = V 2 ahuj/(A + a 2 )^. (46) 
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Our method gives binding energies in units of Huj for sin- 
gle polarons in the constant bare-mass model of 0.2\V 2 a, 
0.32V 2 a and 0.40 V 2 a for a=4, 2 and 1, which may be 
compared with 0.22V 2 a, 0.35V 2 a and 0.45V" 2 a for the 
Holstein model for weak coupling from Eq. (|46|) . Thus 
the fractional changes in binding energies at weak cou- 
pling due to use of the Holstein model appear to be ap- 
proximately 0.05, 0.09 and 0.12 for the three values of a 
considered. The results presented above are for single- 
polaron theory, and what changes would be expected for 
bipolarons is not obvious. However, if we take the results 
based on single-polaron theory seriously for bipolarons, 
then we expect, e.g. for a = 4, V 2 = 0.125, that the 
magnitude of the difference in the energy of the bottom 
of the bare and bipolaron bands to be increased from 
about Q.29Hu> for the constant bare-mass model to be- 
tween 0.30 and O.Slhuj for the Holstein model, compared 
with 0.335huj obtained by El Shawish. Thus in this case 
it appears probable that errors due to the restriction of 
our method to weak or intermediate coupling and due 
to use of a two-parameter trial function for relative mo- 
tion may be of the order of 10% for bipolaron energies at 
Q = 0. 

Besides differences from the Hubbard model, we also 
discuss briefly how far our method is likely to be fairly 
accurate within the framework of the constant-bare-mass 
model that we have used. We note the following points: 

1. From the results mentioned after Eq. (|46|) and values 
of parameters used for the figures and other values men- 
tioned in the first paragraph of this section, the single- 
polaron binding energy lies between about 0.105fio; (for 
V 2 = 0.125 and a = 4) and 0.315^ (for V 2 = 0.375 
and a = 4). The type of variational method we use is 
accurate to about 4% for single-polaron binding energies 
in the Frohlich model and to within 20% for masses^ 
up to polaron couping constants a = 3, corresponding to 
polaron binding energies of about 3huj. Our binding en- 
ergies are so far below this value that our method can be 
expected to be accurate for single polarons, and also for 
bipolarons if one assumes that the coupling constant for 
bipolarons to have the same percentage errors in masses 
as for single polarons is at least half as large as that for 
single polarons. (See point 3 below for further discussion 
of this). 

2. There is also the question of whether we are any- 
where near a transition between large and small polarons. 
Toyozawa^ was the first to study such transitions, for 
the case of interactions between electrons and acousti- 
cal phonons, and found a sudden transition. However, 
since we are dealing with optical phonons here, his work 
does not have much relevance for our problem. For inter- 
actions with optical phonons such transitions (between 
large and nearly-small polarons) were first discussed by 
one of the present authorsi 70 i 71 

Emin 72 considers transitions between large and small 
polarons for a 3-D molecular-crystal model. For the case 
where the bare half bandwidth is ten times as large as 
the phonon energy, he finds a transition between types 



for the lowest-energy state as a function of the binding 
energy of a small polaron which would exist in the case 
of zero bare bandwidth. He calls this energy E b . From 
his Fig. 2 one finds that the transition at T = for the 
above value of bare half-bandwidth occurs when the po- 
laron binding energy (not E b ) is about 2fkv, although 
both types of solutions exist for weak-coupling polaron 
binding energies between about fkv and 3fuv. Even the 
lower bound on these values is considerably larger than 
the values of polaron binding energies of up to 0.315hu> 
which we have discussed in the present paper. Thus it is 
probable that our results are not influenced by any prox- 
imity to such a transition. However, Emin does not dis- 
cuss in detail results in the opposite limit which we have 
considered when bare band half bandwidths are small 
compared with phonon energies. 

For a model with Frohlich electron-phonon interac- 
tions, it appears that a large to small polaron transition 
occurs for coupling constants a near 3 to 5 > 73 i 74 with de- 
tails depending on the degree of adiabaticity. The transi- 
tion is fairly sharp for bare half bandwidths greater than 
phonon energies, but more gradual for the opposite case. 
For the smallest bare half bandwidth of hio considered in 
Refill, their Fig. 3 indicates only a fairly small departure 
of masses from a linear dependence on coupling constant 
a up to one. The cases we concentrate on are a = 4 and 
a = 2, corresponding to bare half bandwidths of 0.31?ia; 
and 1.23huj for a constant bare mass up to the edge of 
the Brillouin zone at it/a, or to 0.125Su and 0.5^ if we 
convert masses to obtain transfer integrals t for a tight- 
binding model with the same bare mass at the bottom 
of the band. Thus, although our coupling constants cor- 
respond to polaron binding energies considerably smaller 
than hw in the Frohlich model, we could not completely 
rule out a larger mass rise than we find by our method 
because of the beginning of a transition between polaron 
types if our band structure were similar to that of a tight- 
binding model as considered in Refs 17311741 However, in 
our case we are assuming that any effect of the lattice 
periodicity on the bare electrons is small, and so in such 
a case any transitions between polaron types will corre- 
spond to transitions in the continuum model, i.e. to a 
change from the lattice following the instantaneous po- 
sition of the electron for weak coupling to responding 
to some average position for stronger coupling, and the 
effects of such transitions are already included in calcu- 
lations such as those reported in Refj6q. 

For our case where there is no significant effect of pe- 
riodicity of the potential on the bare-electron wave func- 
tion, there is no such thing as a small polaron in the 
sense of a state which is a linear combination of states 
with the electron on one lattice site and surrounded by 
the appropriate lattice polarisation, since lattice sites are 
almost indistinguishable from positions in between them. 
In this case the fact that our single-polaron radius for 
a = 4 is only a quarter of the lattice constant does not 
imply small polarons in the usual sense. The only way 
that the lattice constant comes into our model is by a 
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cut-off in the phonon wave vector. Thus in this sense 
it is similar to the continuum-polarisation model with a 
cut-off considered by Schultz.— 

3. Iadonisi et ali2£ also consider transitions between 
types for bipolarons. For a case shown in Fig.l of 
their paper, corresponding to a bare bandwidth equal 
to twenty times the phonon energy, the transition for 
bipolarons occurs at about an 8% smaller values of the 
coupling constant than for polarons. Thus our guess in 
point 1 above that a given percentage error of mass may 
occur at a value of coupling constant of about half that 
of polarons may be pessimistic. 

We get no confirmation in this work of our conjectures 
based on perturbation theory of great enhancements of 
pair binding energies at certain large centre-of-mass wave 
vectors ! 52 ' 53 Also, unpublished calculations of El Shawish 
up to bipolaron wave vectors of ir/a do not give us much 
reason to expect that suggestions of a cusp-like minimum 
at 2tt /a (in an extended zone scheme) indicated by early 
attempts to extend results of our variational method to 
the Hubbard-Holstein model as reported in Ref. [75|, are 
likely to occur in accurate calculations. However, we still 
cannot rule out the possibility that a dip in bipolaron en- 
ergies would be obtained near certain wave vectors if we 
were to use a different type of variational wave function 
for relative motion which could take better advantage of 
the small denominators in the integrands in Eq. (J27J) for 
suitable values of Q and fairly weak coupling. 

Our calculations indicate that there are parameter val- 
ues where bipolarons do not have excessively high masses 
while having binding energies with respect to two single- 
polaron energies greater than a few tenths of the rele- 
vant boson energy. If the bosons are plasmons of en- 
ergy of the order of 2 eV, then this permits bipolaron 
binding energies of the order of 0.5 eV without too great 
increases in masses, whereas for phonons of energies of 
about 0.36 eV, bipolaron binding energies of 0.1 eV can 
be obtained without too large mass increases. A bind- 
ing energy of at least 0.1 eV is a minimum requirement 
for room-temperature superconductivity, assuming pair 
binding energies must be at least about 4fc B T for super- 
conductivity at temperature T. The masses must not be 



too high in order to be able to have a high Bose-Einstein 
condensation temperature for bipolarons without exces- 
sively high bipolaron concentrations. Previous calcula- 
tions of condensation temperatures for bosons with a 
quadratic E(Q) curved in arrays of nanofilaments have 
recently been extended^ 2 - to cases with a dispersion ap- 
proximated by a sum of linear and quadratic terms, as 
indicated to occur for Cooper pairs^ 3 - We hope to modify 
the calculations of Ref. [62| soon by use of a Bogoliubov- 
type of dispersion for pairs (see e.g. Ref. [73), which we 
now think is more appropriate than that based on a 
Cooper-pair model for the strongly coupled pairs at fairly 
low carrier concentrations in which we were interested in 
Ref. M- 

V. CONCLUSIONS 

No support comes from our variational method for pre- 
vious results based on perturbation theory that great en- 
hancements of binding energies of pairs can be obtained 
at appropriate large centre-of-mass velocities. However, 
parameters have been found such that bipolaron masses 
are smaller than about 3m e while keeping binding en- 
ergies with respect to energies of two single polarons 
greater than 0.1 eV. Thus bipolarons in one dimen- 
sion may provide a basis for an explanation of probable 
room-temperature superconductivity in narrow channels 
through films of oxidised atactic polypropylene and other 
polymers, but for different reasons than conjectured in 
earlier papers. 
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